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ABSTRACT 

The clustering properties of local, Si.4GHz > 1 rnJy, radio sources are investigated for 
a sample of 820 objects drawn from the joint use of the FIRST and 2dF Galaxy Red- 
shift surveys. To this aim, we present 271 new 6j < 19.45 spectroscopic counterparts 
of FIRST radio sources to be added to those already introduced in Magliocchetti et 
al. (2002). The two-point correlation function for the local radio population is found 
to be entirely consistent with estimates obtained for the whole sample of 2dFGRS 
galaxies. From measurements of the redshift-space correlation function £(s) we derive 
a redshift-space clustering length s = 10. Ttgi? Mpc, while from the projected correla- 
tion function S(r T ) we estimate the parameters of the real-space correlation function 
£(r) = (r/r y~<, r = 6.7+H Mpc and 7 = 1.6±0.1, where h = 0.7 is assumed. Differ- 
ent results are instead obtained if we only consider sources that present signatures of 
AGN activity in their spectra. These objects are shown to be very strongly correlated, 
with ro = 10.9^2 Mpc and 7 = 2 ± 0.1, a steeper slope than has been claimed in 
other recent works. No difference is found in the clustering properties of radio- AGNs 
of different radio luminosity. Comparisons with models for £(r) show that AGN-fucllcd 
sources reside in dark matter halos more massive than ~ 10 13 4 Mq, higher the cor- 
responding figure for radio-quiet QSOs. This value can be converted into a minimum 
black hole mass associated with radio-loud, AGN-fuelled objects of Mg^ ~ 10 9 Mq. 
The above results then suggest - at least for relatively faint radio objects - the ex- 
istence of a threshold black hole mass associated with the onset of significant radio 
activity such as that of radio-loud AGNs; however, once the activity is triggered, there 
appears to be no evidence for a connection between black hole mass and level of radio 
output. 
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1 INTRODUCTION 

The last twenty years have shown radio sources to be ex- 
tremely good probes of cosmological large-scale structure 
up to significant (z ~ 4) redshifts. Clustering in radio cat- 
alogues was detected in early wide-area studies: Seldner & 
Peebles (1981) and Shaver & Pierre (1989) reported the de- 
tection of slight clustering of nearby radio sources, while 
Kooiman, Burns & Klypin (1995) and Loan, Wall & La- 
hav (1997) measured strong anisotropy in the distribution 
of bright radio objects from the 4.85 GHz Green Bank 
and Parkes-MIT-NRAO surveys. However, it was only few 
years ago that the latest generation of radio surveys such as 
FIRST (Faint Images of the Radio Sky at 20 cm; Becker, 
White & Helfand 1995), WENSS (Rengelink et al. 1998) 
and NVSS (NRAO VLA Sky Survey; Condon et al. 1998) 
included enough objects to allow for high- precision cluster- 
ing measurements (Cress et al. 1996; Rengelink et al. 1998; 
Magliocchetti et al. 1998; Blake & Wall 2003; Overzier et al. 
2003). All these recent analyses reveal the tendency for radio 
objects to be more strongly clustered than optically-selected 
galaxies. 

Assessing the real clustering signal is nevertheless not 
an easy task when it comes to radio sources, since the lack of 
known redshifts for the overwhelming majority of these ob- 
jects forces one to deal with projected quantities such as the 
angular two-point correlation function w(9); unfortunately, 
the relation between angular and spatial measurements is 
dependent on the radio source redshift distribution N(z) - 
which is highly uncertain at the mjy flux levels probed by 
these new surveys (see e.g. Dunlop & Peacock 1990; Maglioc- 
chetti et al. 1999) - and on the unknown redshift evolution 
of the radio clustering signal. This leads to estimates for the 
comoving correlation length ro (here we are assuming the 
spatial two-point correlation function to be described by the 
power-law £(r) = (?"/ro)~ 7 ) which may span the relatively 
broad range ~ 7 — 15 Mpc even for angular measurements in 
good agreement with each other (Magliocchetti et al. 1998; 
Blake & Wall 2003; Overzier et al. 2003). 

Since the natural solution to the above problem (red- 
shift acquisition for all the objects included in a chosen cat- 
alogue) is out of reach with current instruments, a good 
starting point towards the understanding of the clustering 
properties of radio sources is the analysis of homogeneous 
subsamples of objects, bright enough in the visible band to 
allow for optical and spectroscopic follow-ups. 

A first attempt in this direction was performed by Pea- 
cock & Nicholson (1991) who measured the redshift-space 
correlation function for a sample of 310 radio galaxies with 
z < 0.1 and radio fluxes S > 0.5 Jy at 1.4 GHz. These 
authors indeed found radio sources in their catalogue to be 
strongly clustered, with a redshift-space correlation length 
s ^ lift -1 Mpc. 

The present paper analyses the clustering properties of 
~ 820, z < 0.3, Si.4GHz > 1 mjy radio galaxies drawn from 
the joint use of the FIRST and 2dF Galaxy Redshift surveys 
as illustrated in Magliocchetti et al. (2002). By doing this, we 
not only extend the Peacock & Nicholson (1991) measure- 
ments to a statistically more significant sample involving 
less local objects, but we also probe much lower flux den- 
sities where the population contains radio-emitting sources 



that differ from typical AGNs (such as galaxies undergoing 
intense star formation). 

In addition to the above analysis, we will also estimate 
the two-point correlation function (both in redshift space 
and real space) for the homogeneous sample of radio AGNs 
(sources that present signatures of AGN activity in their op- 
tical spectra). This will enable us to use up-to-date models, 
which connect galaxy formation properties with clustering 
behaviour in order to investigate the nature of the dark mat- 
ter halos in which these sources reside. By means of the same 
approach we can also obtain precious information on the na- 
ture of the black hole associated with a radio-active AGN 
to be compared with conclusions from other works. 

For instance, Auriemma et al. (1977) found that radio 
emission is more common in brighter/more massive ellip- 
ticals. This result was recently confirmed e.g. by findings 
of Falomo et al. (2003) , while authors such as Laor (2000) , 
Lacy et al. (2001), McLure & Dunlop (2002), McLure & 
Jarvis (2002) and Dunlop et al. (2003) have investigated the 
connection between black hole mass and radio luminosity, 
coming to conclusions in disagreement with each other. 

The layout of the paper is as follows. Section 2 in- 
troduces 271 new spectroscopic counterparts for local, 
5*1.4 GHz > 1 mjy, radio sources, derived from the joint use 
of the FIRST and 2dF Galaxy Redshift (2dFGRS) surveys. 
Section 3.1 describes some of the properties of the total 
FIRST-2dFGRS sample, which puts together objects pre- 
sented in Section 2 with those introduced in Magliocchetti 
et al. (2002), while Sections 3.2 and 3.3 are respectively de- 
voted to the analysis of the redshift-space correlation func- 
tion £(s) and of the projected correlation function H(r T ) 
both for the whole radio-2dFGRS sample and only for the 
536 radio- AGNs. Section 4 compares the available data with 
models in order to derive the minimum dark matter mass of 
a halo able to host an AGN-fuelled radio source and gives 
some constraints on the minimum black hole mass of a radio- 
loud AGN. Finally, Section 5 summarizes our conclusions. 

Unless stated otherwise, throughout this work we will 
assume Qo = 0.3, A = 0.7, h = 0.7 and as = 0.8 (with as the 
rms density fluctuation within a sphere with a radius of 8ft _1 
Mpc), as the latest results from the joint analysis of CMB 
and 2dFGRS data seem to indicate (see e.g. Lahav et al., 
2002; Spergel et al. 2003). To avoid confusion, we especially 
emphasise that distances will thus normally be quoted in 
Mpc, assuming h = 0.7. Units of h~ x Mpc will occur only 
rarely. This assumption is also made in all figures quoted for 
radio luminosities and black-hole masses. All the correlation 
lengths are in comoving units. 



2 MORE SPECTROSCOPIC COUNTERPARTS 
FOR FIRST RADIO SOURCES 

The 2dFGRS (Colless et al. 2001, 2003) is a large-scale sur- 
vey selected in the photometric bj band from the APM cat- 
alogue (Maddox et al. 1990a, 1990b, 1996) and its subse- 
quent extensions (Maddox et al., in preparation). The final 
version includes 221,414 unique, reliable (i.e. with quality 
flags Q > 3, Colless et al. 2001, 2003) galaxy redshifts to 
an extinction-corrected limit for completeness of bj ~ 19.45 
(with small variations in magnitude limit as a function of 
position over the sky) over an area of 2151 square de- 
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Figure 1. Projected distribution of radio sources with optical counterparts in the APM catalogue for 6j < 19.45. Filled circles identify 
objects with spectral identifications, while empty ones are for those not included in the 2dFGRS. 



1000 



100 



10 



T 



\ r=1000 N - 



s_ o 
° °6 so • 



o 
o 




1000 



100 



10 



r=1000 



10 



_J I I L 



J I "l L 



14 



16 



18 



20 



Figure 2. 6j magnitudes versus radio flux S at 1.4 GHz for different classes of objects. The left-hand panel shows the case for early-type 
galaxies (represented by empty circles), E+AGN (filled circles) and Seyfert galaxies (crosses), while the right-hand panel is obtained for 
the populations of late- type galaxies (filled triangles) and starbursts (filled squares). Dashed lines correspond to constant values of the 
radio-to-optical ratio r = 10, 100, 10 3 , 10 4 ; the dotted lines indicate the threshold of r ~ 30 for radio-loudness. 



grees. The survey geometry consists of two broad declination 
strips, a larger one in the SGP covering the area 3 h 30 m < 
RA(2000) < 21 h 40 m , -37.5° < dec(2000) < -22.5° and 
a smaller one set in the NGP with 9 h 50 m < RA(2000) < 
U h 50 m , 2.5° < dec(2000) < -7.5°, plus 100 random 2- 
degree fields spread uniformly over the 7000 square degrees 
of the APM catalogue in the southern Galactic hemisphere. 
The median redshift of the galaxies is 0.11 and the great 
majority have z < 0.3. 

The completion of the 2dFGRS has allowed us to obtain 
271 new spectroscopic counterparts for Si.4GHz > 1 mjy ra- 



dio objects to be added to the sample presented in Maglioc- 
chetti et al. (2002). As in Magliocchetti et al. (2002), the 
parent radio dataset comes from matching together sources 
in the FIRST (Becker et al. 1995) and APM catalogues over 
the region of the sky between 9 h 48 m < RA(2000) < 14 ,l 32 m 
and -2.77° < dec(2000) < 2.25°, where these two surveys 
overlap. Magliocchetti & Maddox (2002) find 4075 identifi- 
cations - corresponding to 16.7 per cent of the original radio 
sample - in the APM catalogue for 6j < 22 and for a match- 
ing radius of 2 arcsec, value chosen as the best compromise 
to maximize the number of real associations (expected to be 
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~ 97% of the FIRST- APM dataset), while limiting the num- 
ber of random coincidences to a negligible ~ 5 per cent. 971 
objects (hereafter indicated as the photometric catalogue) 
in the parent dataset exhibit b } magnitudes brighter than 
19.45, the limit of the 2dFGRS. 

We note that, even though the fine angular resolution of 
the FIRST survey implies that some of the flux coming from 
extended sources could be either resolved out or split into 
two or more components, leading to a systematic underesti- 
mate of the real flux densities of such sources, this effect has 
been partially corrected for by using the method developed 
by Magliocchetti et al. (1998) to combine multi-component 
objects. The above technique, together with the matching 
procedure, associates optical counterparts to the centroids 
of multi-component sources, thus minimizing the chances of 
missing an object with extended radio emission. Also, as al- 
ready discussed in Magliocchetti et al. (2002) , the correction 
to radio fluxes as measured by FIRST is ~ 1 mjy in the case 
of sources brighter than 3 mjy, and about 30 per cent for 
sources with lmjy < 51.4GHz < 3 mjy (which mainly cor- 
respond to compact spiral galaxies and starbursts). These 
corrections have been shown not to affect any of the results 
obtained in both the previous (Magliocchetti et al. 2002) 
and current works. 

The new spectroscopic counterparts were identified by 
searching in the 2dFGRS catalogue for objects not yet 
present in the Magliocchetti et al. (2002) sample with po- 
sitions that differed by less than 2 arcsec (value of the di- 
ameter of each 2dF fibre) from positions of sources in the 
parent photometric dataset. The properties of these 271 ob- 
jects are described in Table 2 at the end of the paper; for 
each of them the Table indicates: 
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Figure 3. Distribution of monochromatic radio power at 1.4 GHz 
for those FIRST radio sources identified in the 2dFGRS survey. 
The top panel only includes objects that present signatures of 
AGN activity in their optical spectra, while the bottom panel is 
for the whole FIRST-2dFGRS sample. 



(1) Source number 

(2) Right Ascension a (J2000) and (3) Declination 5 (J2000) 
as measured in the 2dFGRS. 

(4) Offset (expressed in arcsecs) between radio and optical 
counterpart in the APM catalogue. 

(5) Radio-flux density (in mjy units) at 1.4 GHz. 

(6) Apparent bj and, when present, (7) R magnitudes of the 
optical counterpart. 

(8) Redshift. 

(9) Spectral Classification. 

Following Magliocchetti et al. (2002), classes for the op- 
tical counterparts of radio sources (column 9 of Table 2) have 
been assigned on the basis of their 2dF spectra. Spectra 
have been compared with known templates (see e.g. Kenni- 
cut 1992; McQuade et al. 1995) which allowed galaxies to 
be divided into 6 broad categories: 

(1) Early-type galaxies, 150 sources where spectra were 
dominated by continua much stronger than the intensity of 
any emission line. These objects can be further divided into 
two sub-classes: 

(1) galaxies with absorption lines only. 

(ii) galaxies with absorption lines + weak [Oil] and Ha 
emission lines denoting little star-formation activity. 

(2) E+AGN-type galaxies, 40 sources showing spectra 
typical of early-types plus the presence of (narrow) emission 
lines such as [Oil], [OIII], [Nil] and [SII], which are strong 
if compared to any Balmer line in emission and indicate the 



presence of large, partially ionized transition regions as is 
the case in active galaxies. 

(3) Late-type galaxies, 51 sources where spectra show 
strong emission (mainly Balmer) lines characteristic of star- 
formation activity, together with a detectable continuum. 

(4) Starburst galaxies, 15 sources with optical spectra 
characterized by an almost negligible continuum with very 
strong emission lines indicating the presence of intense star- 
formation activity. 

(5) Seyfert 1 galaxies, 1 source with spectrum showing 
strong, broad emission lines. 

(6) Seyfert 2 galaxies, 9 sources where the continuum is 
missing and spectra only show strong narrow emission lines 
due to the presence of an active galactic nucleus. 

Distinctions between different classes of sources and in 
particular between E+AGN, Seyfert 2 and Late- type galax- 
ies have relied on the diagnostic emission line ratios of 
Veilleux & Osterbrock (1987), Woltjer (1990) and Rola Ter- 
levich & Terlevich (1997). Note that a definite classification 
was not possible for all the cases. This simply reflects the fact 
that it is in general quite common to find 'composite' galax- 
ies containing both an AGN and ongoing star formation (see 
e.g. Hill et al. 2001). We also found two intermediate cases 
of Seyfert 2-spectral-types with mildly broadened emission 
lines. 

More details on the properties of the full FIRST- 
2dFGRS sample are given in Section 3.1. 
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Figure 4. Left-hand panel: redshift distribution N(z) for the whole spectroscopic sample (solid line) and for the normalized random 
catalogues (dashed line). Right-hand panel: same as before but for the class of radio- AGN sources (see text for details). 



3 THE CLUSTERING PROPERTIES OF 
LOCAL RADIO SOURCES 

3.1 The Data Set 

The combination of the sample introduced in Magliocchetti 
et al. (2002) and of those sources presented in Section 2 
provides us with a total of 828 FIRST radio objects with 
a spectroscopic counterpart in the 2dFGRS. This number 
is reduced to 820 if we only consider spectra with quality 
flag Q>3 and discard the (four) stars present in the dataset. 
Their projected distribution onto the sky is shown in Fig. 1, 
whereby empty circles identify the 971 bj < 19.45 and S > 
1 mjy sources drawn from the parent photometric catalogue 
(see Section 2), while filled dots represent the 820 objects 
with available redshift estimates and spectral classifications. 

In order to investigate the radio properties of the sam- 
ple, Figs 2 and 3 respectively show bj magnitudes versus ra- 
dio flux S at 1.4 GHz and radio power distribution. In more 
detail, the left-hand panel of Fig. 2 is devoted to the b 3 — S 
distribution of "classical" , AGN-powered sources (hereafter 
called radio- AGNs), while the right-hand panel was obtained 
for the sub-class of star-forming galaxies which owe their 
radio emission to processes different from accretion onto 
a central black hole. Empty circles then are for early-type 
galaxies, filled dots for E+AGN, crosses for Seyfert galax- 
ies (regardless of whether type 1 or 2), while filled triangles 
are for late-type galaxies and filled squares for starbursts, 
where the classification follows the one introduced in Sec- 
tion 2. The dashed lines in both panels indicate the loci of 
constant radio-to-optical ratios r = S x lO^J -12 - 5 )/ 2 - 5 ^ and 
the dotted lines illustrate the threshold of r ~ 30 for an ob- 
ject to be considered as radio-loud. It is clear from the plots 
that these two classes of objects tend to occupy different 
regions of the bj — S plane, star-forming sources being in 
general brighter in magnitude and fainter in radio flux. In 
spite of this difference, we see that the overwhelming major- 
ity of the spectroscopic sample has r > 30, i.e. the sample is 



made of radio-loud sources, regardless of the nature of the 
radio emission. 

Radio luminosities have then been derived according to 
the relation V = Pi.4GHz = <Si.4GHz-D 2 (l + z) i+a , and are 
expressed in [W Hz" 1 sr _1 ] units. In the above formula, D is 
the angular diameter distance and a is the spectral index of 
the radio emission (S(u) oc v~ a ). Following Magliocchetti 
et al. (2002), we assumed a = 0.5 for Seyfert 1 galaxies, 
a = 0.75 for early-type galaxies (with or without emission 
lines due to AGN activity), a = 0.7 for Seyfert 2's and 
a — 0.35 both for late-type galaxies and starbursts. Fig. 
3 shows the resulting distribution of radio luminosities for 
the whole spectroscopic sample (lower panel) and for the 
sub-class of radio- AGNs (top panel). The lack of objects at 
logio'P < 21 — 21.5 seen in the distribution of AGN-fuelled 
sources is not seen in the whole spectroscopic sample because 
of the population of low-luminosity star-forming galaxies 
which are limited to radio powers logio'P 22 (see also 
Magliocchetti et al. 2002). We also note that all the sources 
in our sample have moderate luminosities (P < 10 24 ), i.e. 
that all the AGN-fuelled objects belong to the class of FRI 
galaxies (Fanaroff & Riley 1974). This is as expected, since 
the local radio luminosity function of steep spectrum sources 
drops rapidly at V > 10 25 (Dunlop & Peacock 1990), which 
is the typical minimum luminosity of an FRII galaxy. 

With the aim of studying the clustering properties of the 
population of local radio sources, following Hawkins et al. 
(2003) and Madgwick et al. (2003) we have discarded from 
the original FIRST-2dFGRS sample all the radio objects 
in 2dFGRS fields with < 70 % completeness. Furthermore, 
we have also only considered those sources with redshifts 
0.01 < 2 < 0.3 and magnitudes b 3 < 19.37. This magnitude 
limit is the brightest over the area of the 2dFGRS survey 
(Colless et al. 2001, 2003), and so using this limit over the 
whole area removes the need to correct for variations in the 
magnitude limit. These cuts leave us with a total of 761 
objects (spectroscopic sample) out of which 536 belong to 
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the population of AGN-fuelled sources (radio- AGN sample) . 
Their redshift distributions are illustrated by the histograms 
in Fig. 4. 

3.2 The Redshift-Space Correlation Function 

The standard way to quantify the clustering properties of 
a particular class of sources is by means of the two-point 
correlation function £ which measures the excess probabil- 
ity of finding a pair in the two volume elements dV\ and 
dV2 separated by a distance r. In practice, £ is obtained by 
comparing the actual source distribution with a catalogue of 
randomly distributed objects, subject to the same redshift 
and mask constraints as the real data. We chose to use the 
estimator (Hamilton 1993) 



as) = 4 x 



DP RR 
{DRY 



1, 



(1) 



where DD, RR and DR are the number of data-data, 
random-random and data-random pairs separated by a dis- 
tance s (note that throughout this section we will write the 
distance as s instead of r to stress the fact that we are work- 
ing in redshift space). 

We generated random catalogues with ten times as 
many objects as the real datasets, and modulated the an- 
gular distribution by both the 2dFGRS and FIRST cover- 
age maps, so that the instrumental window functions did 
not affect the measured clustering. The redshift distribu- 
tion for these objects was then generated using the se- 
lection function of the different datasets under considera- 
tion. Equivalent Wi(z) = 1/[1 + 47r tir(z)J3(s)] = 1 (with 
ur(z) number density of radio sources at redshift z and 
J 3 = /J £(s')s' 2 ds'; see Hawkins et al. 2003 and Madg- 
wick et al. 2003) weights were assigned to each object 
both in the real and random catalogues; this assumption 
was justified by the extremely low number density of ra- 
dio sources (n R (z) max ~ 6 10~ 6 Mpc~ 3 ) which makes the 
47r nn{z)J:i{s) contribution to w» a negligible one for any 
sensible value of J3. 

The above procedure was repeated both for the whole 
spectroscopic sample and in the case of only radio-AGNs. 
The redshift distributions for the two random catalogues 
corresponding to the different samples are shown by the 
dashed lines in Fig. 4. The resulting redshift-space correla- 
tion functions £(s) are shown in Fig. 5, where the left-hand 
panel is for the full spectroscopic sample, and the right-hand 
panel for radio-AGNs. The error bars on £(s) were obtained 
by jack-knife resampling each data-set 20 times. 

Fig. 5 suggests that a universal power-law is not a good 
fit on all scales for either of the measured £(s)'s and that it 
is particularly poor for the AGN-powered sources. Although 
a power-law fits well over the range 6 < s/[Mpc] < 50, the 
redshift-space correlation function steepens at larger scales, 
and flattens on smaller scales, particularly for the radio- 
AGN sample. We estimate the true value of the redshift- 
space correlation length so by fitting a localized power-law of 
the form £(s) = (s/so)~ 7s only in the range 6 < s/[Mpc] < 
50. 

The corresponding best values for the parameters so 
and 7s can be obtained by a x fit to £ as a function of s. 
This analysis gives s = 10.7jlo ® M P C and 7= = 1.5 ± 0.1 
for the whole FIRST-2dFGRS sample, while in the case of 



AGN sources one gets s = 13 ± 0.9 Mpc and 7 S = l.Sj:^. 
These best-fits are shown by the dashed lines in Fig. 5. 
The dotted line on the left-hand panel of Fig. 5 shows the 
redshift-space correlation function for the entire population 
of 2dFGRS galaxies, derived by Hawkins et al. (2003); the 
similarity between this measurement and the one for local 
radio-sources is remarkable. This can be understood because 
the FIRST-2dFGRS catalogue probes radio fluxes down to 
Siaghz = 1 mjy, which results in a mixture of "classical" 
AGN-fuelled sources, generally hosted by elliptical galaxies, 
and star-forming objects such as spirals and irregulars (see 
Section 2). This mix of early and late type galaxies makes 
the dataset look similar to a fair (even though very sparse) 
sample of the whole population of 2dFGRS galaxies. 

The situation however looks very different if one consid- 
ers only AGN-fuelled sources. In this case the sample is much 
more strongly correlated; the correlation length is higher, 
and also £(s) has a steeper slope. We will discuss the im- 
plications of these findings in Section 3.3 when dealing with 
the real-space correlation function. 

As a final step, we are interested in investigating 
whether there is any dependence of the clustering signal on 
radio luminosity. We divided the radio- AGN sample into two 
distinct subsets depending on whether the sources have ra- 
dio luminosity brighter or fainter than logioP = 22.0. There 
are 260 objects in the faint sample, and their average radio 
luminosity is logioP = 22.9. There are 276 objects in the 
bright sample, with average radio luminosity logioP = 21.6. 

The redshift-space correlation functions for the two 
samples are plotted in Fig. 6; the fainter sample is plotted as 
filled circles, and the brighter sample as open squares. Even 
though the small number of sources in the subsets lead to 
large uncertainties in the measurements, it is clear that the 
two £(s)'s are entirely compatible with each other. The best- 
fit parameters are so = 12 ± 3 Mpc and 7 S = 1.5 ± 0.15 for 
brighter AGNs and s = 13 ± 3 Mpc and 7 S = 1.4 ± 0.1 for 
the fainter ones. This result shows no evidence for luminos- 
ity dependence in the clustering amplitude, and shows that 
the amplitude of the correlation function at 10 Mpc cannot 
differ by more much than a factor of 2 even though the two 
samples differ by a factor of 20 in luminosity. This is con- 
sistent with the results of Peacock & Nicholson (1991) who 
analysed an all-sky sample of radio galaxies at z < 0.1, and 
found no evidence for any difference in clustering amplitude 
for samples 22.5 < logio"P < 23.5 and 23.5 < logioP < 24.5. 
These results imply that the correlation amplitude of radio 
galaxies is independent of their radio luminosity, at least in 
the power range 10 20 < V < 10 25 . 

3.3 The Projected Correlation Function 

The previous Section has shown that the redshift-space 
correlation function is extremely useful for characterizing 
the main properties of the clustering of low-redshift radio 
sources. However, £(s) differs from the real-space correlation 
function because of peculiar velocities that lead to redshift- 
space distortions (see e.g. Hawkins et al. 2003 and Madgwick 
et al. 2003). It follows that both the amplitude and the slope 
of the real-space correlation function f (r) will differ from the 
values measured in Section 3.2. 

The standard way to recover £(r) is to first compute 
the two-dimensional correlations as a function of separation 
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Figure 5. Left-hand panel: redshift-space correlation function for the whole spectroscopic sample (761 sources). Error bars are obtained 
via jack-knife resampling of the dataset under consideration; given the large uncertainties associated with the smallest-scale measurement 
in the corresponding bin we only plot its upper limit. The dashed line represents the best fit to the data (derived only for those points 
included in the range 6 < s/[Mpc] < 50) corresponding to 7„ = 1.5 and sq = 10.7 Mpc, while the dotted line illustrates the best-fit results 
found by Hawkins et al. (2003) for the entire population of 2dFGRS galaxies. Right-hand panel: redshift-space correlation function for 
the sub-sample of radio- AGNs (536 sources). Error bars are obtained as for the spectroscopic sample. The dashed line represents the 
best fit to the 6 < s/[Mpc] < 50 data, corresponding to 7 S = 1.8 and so = 13 Mpc. 
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Figure 6. Redshift-space correlation function for the two sub- 
samples of radio-AGNs respectively brighter (open squares) and 
fainter (filled circles) than log^T 3 = 22. Error bars are from Pois- 
son noise estimates for the dataset under consideration. 



parallel and transverse to the line of sight, £(r- P ,r T ) and 
then integrate it along the r P direction in order to obtain 
the projected correlation function 3 as: 



(r T ) = 2 



(2) 



We estimate £(r P , r T ) by means of equation (1) - where this 
time pairs are counted in a grid of bins Ar P , Ar T - and 
derive 3 by means of equation (2). We set the upper limit 
for the integral to r P lax = 60 Mpc as the best compromise 
between a large enough value to produce stable results and 
a small enough value to avoid extra noise added to 3. 

As in the previous subsection, 3 has been evaluated 
both using the whole spectroscopic sample and the sample 
of objects showing signatures of AGN in their spectra. Once 
again, the errors were calculated from 20 jack-knife resam- 
plings of the datasets under consideration. 

The resulting measurements on scales 1-30 Mpc are 
shown by the solid points in the two panels of Fig. (7): the 
left panel shows the results for whole spectroscopic sam- 
ple, and right panel show the results for AGN sample. The 
open squares in Fig. (7) show the projected correlation func- 
tions from corresponding optical 2dFGRS samples: in the 
left panel they show the results for the entire population of 
2dF galaxies (Hawkins et al. 2003); in the right panel they 
show the results for early-type 2dFGRS galaxies (Madgwick 
et al. 2003). From these comparisons we see that the cluster- 
ing properties of the whole spectroscopic sample are remark- 
ably similar to those of "normal" galaxies, but the projected 
correlation function of radio AGNs has an amplitude that is 
roughly twice as large as that of local ellipticals. 

From the measurements presented in Fig. 7 it is then 
possible to estimate the real-space correlation function £(r) 
via (Davis & Peebles, 1983) 



3(tt) = 2 / £(r) 



r dr 



(r 2 



r 2\l/2 



(3) 



which, if we assume the power-law form £(r) = (r/ro) 
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Figure 7. Left-hand panel: projected correlation function for the whole spectroscopic sample (761 sources). Error bars are obtained as 
jack-knife resampling of the dataset under consideration. The dashed line represents the best fit to the data corresponding to 7 = 1.6 
and ro = 6.7 Mpc, while open squares illustrate the results obtained by Hawkins et al. (2003) for the entire population of 2dFGRS 
galaxies. Right-hand panel: projected correlation function for the sub-sample of radio-AGNs (536 sources). Error bars are obtained as 
for the spectroscopic sample. The dashed line represents the best fit to the data corresponding to 7 = 2.0 and ro = 10.9 Mpc, while open 
squares illustrate the results obtained by Madgwick et al. (2003) for the population of early-type 2dFGRS galaxies. 



can be integrated analytically leading to 

H(rr) = ry^H,, (4) 



Table 1. Best-fit parameters to £. In the case of £(s) = (s/so) u 
the fit only uses points with 6 < s/[Mpc] < 50 Mpc, while for ro 
and 7 all the points between r T = 1 and r T = 30 Mpc have been 
considered. 



where s T =r(i)r(^i)/r(i). 

A least-square fit to the data for the whole spectroscopic 
sample gives ro = 6.7lj'j Mpc, 7 = 1.6 ± 0.1, surprisingly 
close to the results of Hawkins et al. 2003. For the AGN 
population we find r = 10.9t^ Mpc, 7 = 2.0 ± 0.1, which 
is similar to that found by Norberg et al. (2002) for the 
brightest early- type galaxies in the 2dFGRS. A summary of 
the results derived here and in Section 3.2 can be found in 
Table 1. 

We note that a measured correlation length ro — 
11 Mpc is in very good agreement with previous estimates 
derived from deprojection of the angular two-point corre- 
lation function w(8). These include the results of Cress et 
al. (1996), Loan et al. (1997), Magliocchetti et al. (1998) 
and Magliocchetti & Maddox (2002). However, more recent 
works (e.g. Blake & Wall 2002; Overzier et al. 2003) seem to 
find lower values for ro, of the order of 5-6 h~ Mpc. Given 
that the functional form for the redshift distribution of ra- 
dio sources at the mjy level adopted by the various authors 
is the same, i.e. the one proposed by Dunlop & Peacock 
(1990), we ascribe the differences found by Blake & Wall 
(2002) and Overzier et al. (2003) as caused by not properly 
dealing with the issue of multi-component sources. Ignor- 
ing the effect of multiple-sources on the observed clustering 
properties of radio objects in fact leads (amongst others) to 
an underestimate of the slope of w(9), 1 — 7, which in turn 
yelds smaller inferred values for the corelation length ro. 
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4 CONSTRAINTS ON THE HALO AND 
BLACK HOLE MASS 

Since the H measurements presented in the right-hand panel 
of Fig. 7 have been obtained for a homogeneous sample of 
AGN-fuelled sources, we can easily predict the clustering for 
a range of models to be directly compared with the data. 

Under the assumption that each dark matter halo can 
only host one radio- AGN (but see later for a more general 
discussion), the predicted spatial correlation function £ (on 
a scale r and at a redshift z) for this class of sources can in 
fact be written as: 

i(r,z) = £dm(r,z) 6^ ff (M min ,2). (5) 

This expression holds on scales > 1 Mpc where the halo- 
halo exclusion effects are negligible (for the general case see 
Magliocchetti & Porciani 2003). 

The mass-mass correlation function £dm(r, z) in Equa- 
tion (5) is fully specified for a given cosmological model and 
a chosen normalization as- The precise form of £dm(r, 2) can 
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be analytically derived following the approach of Peacock 
and Dodds (1996). The effective bias & e ff (M m i n , z), which 
determines the way radio sources trace the mass distribu- 
tion, is obtained via: 



&eff(M min ,z) = 



b(m, z) n(m, z) dm 
n(m, z) dm 



(6) 



where b(m, z) and n(m, z) respectively are the linear bias 
factor and the halo mass function of individual halos of mass 
m at a redshift z and - in our specific case - M m i n is the 
minimum mass of a halo able to host a radio- AGN. 

In order to compare models with data we have then 
written both b(m, z) and n(m, z) according to the Sheth 
& Tormen (1999) prescriptions and derived the theoretical 
projected E by integrating the spatial correlation function 
(5) along the direction parallel to the line of sight r P via 
equation (3), where all the relevant quantities are evaluated 
at the median redshift of the 2dFGRS, (z) ~ 0.1. 

The dotted lines in Fig. (8) show the predictions ob- 
tained from equation (5) for different values of the mini- 
mum halo mass necessary to host a radio- AGN as compared 
to the measurements derived in Section 3.3. The chosen val- 
ues for the minimum mass range from M m i n = 10 10 Mq 
(bottom dotted curve) to 10 14 Mq (top dotted curve). De- 
spite the large error bars, it is clear that none of the pre- 
dicted S's for minimum halo masses < 10 13 Mq can pro- 
vide a reasonable fit to the data, since they all fall be- 
low the observed values. Indeed, the best fit is obtained 



13.7 



+0.2 



corresponding to an ef- 



We note that very similar values 



for logio[M min /M ] = 
fective bias & e ff = 1.85. 
(i.e. logio[M min /M ] = 13.5 ±°-j and logic, [M min /M ] = 
13.0lg'g) are also found if one releases the assumption of 
one radio galaxy per halo and considers more general distri- 
butions of the kind (N rad i ) oc m 0,5 (see e.g. Peacock 2003) 
or a more extreme (N ra( jio) oc m (with (N ra dio) average num- 
ber of radio galaxies per halo of mass m) even though, as 
expected, the difference between models derived for different 
values of M m i n decreases as one increases the weight given 
to larger-mass haloes (see Fig. 8). 

The implications of the above results are quite intrigu- 
ing since they suggest, independent of the chosen halo oc- 
cupation form, that radio-AGNs reside in relatively massive 
haloes, spanning a mass range from rich groups to super- 
massive clusters. Smaller halos seem to be inhibited from 
hosting a radio-emitting AGN. This conclusion is also sup- 
ported by the small-scale behaviour of the redshift-space 
correlation function of radio-AGNs as measured in Section 
3.2. Indeed, both the lack of pairs with distance < 1 Mpc 
and the flattening of £(s) on scales < 5 Mpc have a natural 
explanation if one assumes AGN-powered sources to be con- 
fined in massive halos (even if one cannot exclude this latter 
finding as due to the fingers of God effect, see Hawkins et 
al. 2003). In the case where only one source per halo is as- 
sumed, the halo-halo spatial exclusion ensures that no pair 
of radio-AGNs can be found with a separation smaller than 
the sum of the virial radii of the halos that host them. For 
a typical halo mass of 10 13 ' 4 Mq the minimum separation 
for the adopted cosmology turns out to be ~ 1.6 Mpc. Al- 
lowing for the uncertainties in our determination of M m i n , 
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Figure 8. Projected correlation function of radio-AGNs. Mea- 
surements are those obtained in Section 3.3, while the dotted 
lines show the predictions derived from equation (5) for mini- 
mum halo masses ranging from 10 10 (bottom curve) to 10 14 (top 
curve) Mq. Dashed lines are obtained as in the former case but 
for a mean halo occupation number (N radio ) oc m (see text for 
details). 



this implies that there should be a deficit of pairs on scales 
~ 1 — 5 Mpc, consistent with the observed flattening of £(s) 
on small scales. 

Thus our average estimate of the minimum halo mass 
for a radio- AGN is logio [M m i n /MQ] ~ 13.4), which is about 
a factor 10 larger than that associated with "normal" early- 
type galaxies, logio[M m i n /MQ] ~ 12.5, (Magliocchetti & 
Porciani 2003). This is also about a factor 10 larger than 
that obtained for the population of radio-quiet quasars in the 
2dF QSO dataset, logio[M min /M Q ] ~ 12.1 - 12.6 (Grazian 
et al. 2003). All three of the above determinations are based 
on samples drawn from the same parent catalogue (APM 
survey), and so selection effects are not expected to bias the 
comparison between the different results. 
Taking these numbers at face value suggests that radio- 
AGNs and radio-quiet quasars are two distinct populations, 
with more massive dark matter halos required for an AGN 
to trigger radio activity. 

On the other hand, the results presented in Section 
3.2 (see Fig. 6) together with those obtained by Peacock 
& Nicholson (1991) for more luminous radio galaxies show 
that radio-AGNs with very different radio luminosities have 
very similar correlation amplitude. If one thinks in terms 
of the model expressed in equation (5), there appears to be 
little or no correlation between radio luminosity and min- 
imum halo mass. A simple comparison of Figs 6 and 8 il- 
lustrates that the correlation function estimates obtained 
for low-luminosity and moderately bright-luminosity radio- 
AGNs both fall in the same 10 13 < M min /M < 10 14 range 
as found for the whole sample). 

If we then combine the above information, the data 
seem to suggest that there is a threshold halo mass re- 
quired to produce significant radio emission from an AGN 
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but - once the radio activity is triggered - there seems to be 
no connection between radio luminosity of the sources and 
dark matter content of those halos hosting them. Note that 
this statement does not necessarily imply the minimum halo 
mass to be the only factor playing a relevant role in the on- 
set of radio emission from an AGN. Also the threshold may 
not be a sharp mass threshold, but could be as a smooth 
transition between the radio-quiet and radio-loud regimes 
for halos near M m i n . 

These results on the minimum halo mass can also be 
viewed in terms of a minimum black hole mass. Ferrarese 

(2002) gives a simple prescription to convert the dark matter 
mass Mdm of a halo hosting a black hole into black hole mass 
M BH : 

-M™-~0 1(| Mdm (7) 
10 8 M Q \IQ™Mq) { ' 

Thus, a minimum halo mass M m i n ~ 10 13 ' 4 M© corresponds 
to Mj^j ~ 10 9 Mq, which would represent the minimum 
black hole mass required to onset radio-AGN activity - at 
least for relatively faint/FRI sources such as the overwhelm- 
ing majority of those we have considered in our analysis 
(see Section 3.1). It is interesting to note that this estimate 
is in good agreement with both the results of Marchesini, 
Ferrarese & Celotti (2003) and also those of Dunlop et al. 

(2003) , despite the entirely different methods used in our 
analysis. Also our minimum black hole mass is noticeably 
larger than that obtained for radio-quiet quasars; from an 
analysis of the 2dF QSO dataset in the same redshift range 
spanned by the sources considered in this work, Corbett et 
al. (2003) find 6.5 < logio [Mbh/Mq] < 8.3. By making use 
of equation (7), we can then conclude that - at least for 
moderately faint/FRI radio sources - there seems to be a 
minimum black hole mass of about 10 9 Mq required to pro- 
duce significant radio emission from an AGN; however, once 
the radio emission is produced, there is little or no correla- 
tion between radio luminosity and black hole mass. This is 
in agreement with other analyses of radio-loud quasars by 
Dunlop et al. (2003) and Cirasuolo et al. (2003). 

We note that, even though the connection between 
black hole mass and radio power has been obtained in our 
analysis as a consequence of the connection between radio 
power and halo mass, this first relationship is supposed to 
be more fundamental and physically based as it is the black 
hole which actually powers the AGN. Any correlation be- 
tween radio power and halo mass is then expected to arise 
as a secondary effect, since black hole mass is observed to 
scale with bulge luminosity and this is turn is related - in 
the case of elliptical galaxies - to halo mass (see e.g. Magor- 
rian et al. 1998; Kormendy & Gebhardt 2001; McLure & 
Dunlop 2002; Archibald et al. 2002; Granato et al. 2004) 



5 CONCLUSIONS 

This paper has presented the clustering properties of local, 
5i.4GHz > 1 mjy radio sources by making use of a sam- 
ple of 820 objects drawn from the joint use of the FIRST 
survey and the 2dFGRS. To this aim, we have introduced 
271 new bj < 19.45 spectroscopic counterparts for FIRST 
radio sources to be added to those already obtained by 
Magliocchetti et al. (2002). These objects can be divided in 



two broad sub-classes: (1) star-forming objects (66 sources, 
which include both late- type and star-burst galaxies), which 
owe their radio activity to processes connected to intense 
star formation and (2) radio- AGNs (202 sources), where the 
radio signal stems from accretion processes onto a central 
black hole. In both cases, the redshift range range spanned 
extends up to z ~ 0.3 (with star-forming objects being rela- 
tively more local than AGN-fuelled sources) and radio lumi- 
nosities cover the interval 10 20 < V < 10 24 , which identifies 
the radio-AGNs included in the sample as FRI sources. 

The redshift-space correlation function £(s) and the 
projected correlation function H(r T ) have been calculated 
for both the total FIRST-2dFGRS sample and for the sub- 
class of radio-AGNs. The results for the two populations are 
quite different. 

In the case of all FIRST-2dFGRS sources, both the cor- 
relation functions are found to be entirely consistent with 
the estimates obtained by Hawkins et al. (2003) for the 
whole sample of 2dFGRS galaxies. From measurements of 
the redshift-space correlation function £(s) = (s/so)~ 7s we 
derive a redshift-space clustering length so = 10.7io.? Mpc 
and a slope 7 S = 1.5 ± 0.1, while from the projected 
correlation function H(r T ) we estimate the parameters of 
the real-space correlation function £(r) = (r/ro)~ 7 to be 
r = 6.7l ( i;? Mpc and 7 = 1.6 ± 0.1. 

Sources that show signatures of AGN activity in their 
spectra are found to be very strongly clustered, with a value 
for the redshift-space clustering length, so = 13.0 ± 0.9 (in 
excellent agreement with the findings of Peacock & Nichol- 
son 1991). Measurements of the real-space correlation func- 
tion lead to ro ~ 11 Mpc and 7 ~ 2, slightly steeper and 
higher amplitude than optically selected galaxies, but very 
similar to bright early- type galaxies. 

We also find no significant differences in the clustering 
properties of faint (V < 10 22 ) compared to brighter (V > 
10 22 ) AGN-fuelled radio sources. 

Comparisons with physically-motivated models for the 
clustering properties of classes of galaxies show that AGN- 
fuelled sources have to reside in dark matter halos more 
massive than logio [M m ;„/MQ] ~ 13.4, higher than the figure 
recently measured for radio-quiet QSOs (see e.g. Grazian 
et al. 2003). Under certain assumptions, this value can be 
converted into a minimum black hole mass associated with 
radio-loud, AGN-fuelled objects of Mg£ ~ 10 9 M Q , again 
larger than current estimates for the typical black hole mass 
associated with local 2dF (radio-quiet) quasars (Corbett et 
al. 2003). 

The above results then suggest - at least for moder- 
ately faint/FRI radio sources such as those included in our 
sample - the existence of a threshold halo/black hole mass 
associated with the onset of significant radio activity such 
as that of radio-loud AGNs. We stress that such threshold 
is not necessarily a sharp one, as it could as well identify 
an allowed range (in the high-mass regime) for the transi- 
tion between radio-quiet and radio-loud regimes. Once the 
activity is triggered there seems to be no evidence for a con- 
nection between halo/black hole mass and radio luminosity. 
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